Non-existence of stationary two-black-hole configurations 

Gemot Neugebauer ■ Jorg Hennig 



Abstract We resume former discussions of the question, whether the spin-spin repulsion 
and the gravitational attraction of two aligned black holes can balance each other. To answer 
the question we formulate a boundary value problem for two separate (Killing-) horizons 
and apply the inverse (scattering) method to solve it. Making use of results of Manko, Ruiz 
and Sanabria-Gomez and a novel black hole criterion, we prove the non-existence of the 
equilibrium situation in question. 
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1 Introduction 

This paper is meant to contribute to the present discussion about the existence or non- 
existence of stationary equilibrium configurations consisting of separate bodies at rest. Her- 
mann Weyl, whom Jiirgen Ehlers admired especially, was the first person to consider the 
problem of two separate static (axisymmetric) bodies in equilibrium 1271 . To mention only 
one modern advancement in this field we refer to a paper by Beig and Schoen |5 1, who were 
able to prove a non-existence theorem for a reflectionally symmetric static «-body configu- 
ration. 

Our intention is to involve the interaction of the angular momenta of rotating bod- 
ies ("spin-spin interaction") which could generate repulsive effects compensating the om- 
nipresent mass attraction. A characteristic example for such a stationary configuration could 
be the equilibrium between two aligned rotating black holes. We will present and review a 
chain of old and new arguments which finally forbid the equilibrium situation. 
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Fig. 1 Illustration of the two-black-hole equihbrium situation in Weyl coordinates. The event horizons J^i 
and J^2 of the two black holes are located in the intervals [^"2,^^!] and [^Ci.A's] on the ^-axis, respectively. 
The remaining parts iz^'" of the ^-axis correspond to the rotation axis. 



Our argumentation is based on a boundary value problem for two separate (Killing-) 
horizons (see Fig. [TJ and the characterization of sub-extremal black holes by Booth and 
Fairhurst |6| and follows the ideas of Manko and Ruiz 1 19] who solved the equilibrium 
problem for the so-called double-Kerr solution. 

The double-Kerr (more precicely: double-Kerr- NUT) solution, first derived in l lTlllSI , is 
a seven parameter solution constructed by a two-fold Backlund transformation of Minkowski 
space. Since a single Backlund transformation generates the Kerr-NUT solution that con- 
tains, by a special choice of its three parameters, the stationary black hole solution (Kerr 
solution) and since Backlund transformations act as a "nonlinear" superposition principle, 
the double-Kerr-NUT solution was considered to be a good candidate for the solution of 
the two horizon problem and extensively discussed in the literature ll5lll4ll26lll2ll28lfT3l l8l 
fT6lll8lfT7lfT9l . However, there was no argument requiring that this particular solution be the 
only candidate. In this paper we will remove this objection and show that the discussion of 
a boundary value problem for two separate horizons necessarily leads to the double-Kerr- 
NUT solution. Thus we can make use of the equilibrium conditions for this solution which 
ensure that the intervals i/", (see Fig.[T]l are regular parts of the axis of symmetry. 
After a too restrictive ansatz in (15|, Tomimatsu and Kihara 1 14,26] derived and discussed 
a complete set of equilibrium conditions on the axis. Reformulations and numerical studies 
by Hoenselaers |13 ] made plausible that the two gravitational sources (black hole candi- 
dates) of the double-Kerr-NUT solution (located at the intervals Q = 0, Ki > C, > K2 and 

> ^ > ^^4) cannot be in equilibrium if their Komar masses are positive. The first decisive 
step toward prove the Hoenselaers conjecture was taken by Manko, Ruiz and Sanabria- 
Gomez ]18], who derived an explicit and easily applicable form of the Tomimatsu-Kihara 
equilibrium conditions and, as an important complement, analytical formulae for the Komar 
masses and angular momenta of the gravitational sources. Manko and Ruiz fl91 completed 
their non-existence proof by showing that the equilibrium conditions for the double-Kerr- 
NUT solution are indeed violated for positive Komar masses. This is, however, a critical 
point of their analysis. To the best of our knowledge there is no argument in favour of the 
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positiveness of the Komar mass. (On the contrary, Ansorg and Petroff fl] have given a con- 
vincing counterexample.) 

In this paper we replace the Komar mass inequality (positivity of the Komar mass of 
each black hole) by an inequality connecting angular momentum and horizon area |9|. This 
relation is based on the causal structure of trapped surfaces in the interior vicinity of the 
event horizon |6l. In this way we can complete the no-go theorem, avoiding more laborious 
investigations of the domain outside the horizons and off the axis of symmetry (e.g., the 
search for singular rings or other singularities — in Sec. l5.2l we will return to that question). 



2 The boundary value problem 

2.1 The boundary conditions 

The exterior vacuum gravitational field of axisymmetric and stationary gravitational sources 
can be described in cylindrical Weyl-Lewis-Papapetrou coordinates by the line element 

d.?2 = e-^u [(P-''{dg^ + dC^) + g^difi^] - e^^ {dt + ad<p)^ , (1) 

where the "Newtonian" gravitational potential U, the "gravitomagnetic" potential a and the 
"superpotential" k are functions of g and ( alone. Figure [T] shows the boundaries of the 
vacuum region: .e/^, i?/", x/^ are the regular parts of the axis of symmetry, J^i and are 
Killing horizons and 'W stands for spatial infinity. Regularity of the metric along jt/", 
jz/^ means elementary flatness and uniqueness on the axis of symmetry, 

s/^,s/'^: a = 0, k = 0. (2) 

The spacetime has to be fiat at large distances from the horizons, 

U^O, a^O, k^O, (3) 

i.e. the line element (HJ takes a Minkowskian form in cylindrical space (gX,v)- time (t) 
coordinates. 

The metric (|T]| allows an Abelian group of motions G2 with the generators (Killing 
vectors) 

5' = 5;, ^'^i<0 (stationarity), (4) 
ri' = S'^, rj'rii>0 (axisymmetry) , (5) 

where the Kronecker symbols S', and 5^ indicate that ^' has only a ^-component whereas rj' 
points in the azimuthal (/3-direction along closed circles. Obviously, 

e^f^ = -e^i, a = -c-^'^me (6) 

is a coordinate-free representation of the two relativistic gravitational potentials U and a 
with the boundary values ([2li, ll3]l. 

In stationary and axisymmetric spacetimes, the event horizon of a black hole is a Killing 
horizon which can be defined by a linear combination L of ^ and -q, 

L = ^ + nri, (7) 
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where ^7 is a constant. A connected component of the set of points with e^^ := —{L,L) = 0, 
which is a null hypersurface, (de^*^, de^^) = 0, is called a Killing horizon J^{L), 

Jf(L): e2^ = -(L,L)=0, (de2^,de2*') = 0. (8) 

Since the Lie derivative of e^^ vanishes, we have (Ljde^'') = 0. Hence, L and e^^ being 
null vectors on M'{L) are proportional to each other, 

jr(L) : de^^ = -2kL. (9) 

Using the field equations one can show that the surface gravity «: is a constant on M'{L). In 
Weyl-Lewis-Papapetrou coordinates the event horizon degenerates to a "straight line" and 
covers a (^-interval at p = |7 |. To formulate boundary conditions on the horizons M\ and 
(see Fig.[T]l we make use of l|6ll and ^ to express e^'' in terms of e^^, a and g, 

M{: e^^i :=e2^ [(l + «ifl)2-r2?£iV'*^] =0, p = 0, Ki>C>K2, (10) 
Mi- e^^^ ■■=e^" [{l + 02af-nlQ^&''^"]=Q, q = 0, K3>C>K4. (11) 

ill, fi2 are the constant angular velocities of the horizons M{, M2, respectively. 



2.2 The field equations 

The vacuum Einstein equations for the metric potentials U , a, k are equivalent to the Ernst 
equation 

m) [f.ee + f.CC + -Z^) = fl + /,C (12) 

for the complex function 

f{Q,0 = t'''^^'^Hib{e,0, (13) 

where b replaces a via 

fl,e = i?e-'*f^fc^, a^ = -ee-4^^7,g (14) 
and k can be calculated from 

Ke = s[ul-Ul+^-^-'"{b^^^-b\)\, (15) 
fc, C = 2q \u,gU,(^ + ie-4^^.,,Z7,j] . (16) 

As a consequence of the Ernst equation ( 112b . the integrability conditions a^gQ = a^i^g 
and k gQ = k ^g are satisfied such that the metric potentials a and k may be calculated via 
line integration from the Ernst potential /. Thus the boundary value problem for the vacuum 
Einstein equations reduces to a boundary value problem for the Ernst equation. However, we 
have to cope with non-local boundary conditions for the Ernst potential, see l[3]l, l ll4b , 
( I15l >, i ll6b . Fortunately, these boundary conditions are well-adapted to the "inverse method", 
which will be applied to solve the boundary value problem. 
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2.3 The inverse method 



The inverse (scattering) method first apphed to solve initial value problems of special classes 
of non-linear partial differential equations in many areas of physics (Korteweg-de Vries 
equation in hydrodynamics, non-linear Schrodinger equation in non-linear optics etc.) is 
based on the existence of a linear problem (LP) whose integrability condition is equivalent 
to the non-linear differential equation. Luckily, the Ernst equation has an LP too, so one can 
try to tackle boundary value problems for rotating objects including black holes. We use the 
LPHoEa 

r / R n \ 

+ A 



B 
OA 



B 
A 



A 
B 



1 



OA 
B 



(17) 



where #(z,z, A) is a 2 x 2 matrix depending on the spectral parameter 



K- 



K + iz 

as well as on the complex coordinates z = Q + \-C,,z = g — iC,, whereas 

. f,z „ f.z 



f+r 



B = 



f + f 



(18) 



(19) 



and the complex conjugate quantities A, B are functions of z, z (or g, Q alone and do not 
depend on the constant parameter K. From the integrability condition # = # =, and the 
relations 



I] 



1 

4^ 



1) 



(20) 



it follows that the A -independent coefficients of a matrix polynomial in X have to vanish. 
The result is the Ernst equation illi . Vice versa, the matrix # calculated from jlll does 
not depend on the path of integration if / is a solution to the Ernst equation. The idea of 
the inverse (scattering) method is to construct for fixed but arbitrary values of z, z, as 
a holomorphic function of A and to calculate f{g,C) from To obtain the dependence 
on A for the two-horizon problem we have to integrate the linear system ( 117b along the 
dashed line ^ = .c/"UJ^Uii/°Uj^U ssf^ U making use of the boundary conditions 
iHl, (H, (Toll, (HB- As was shown in (23) and IIU (see Eqs. (34), (45), (57), (58) in (lH) the 
result of the integration is a matrix representation of the axis values of the Ernst potential 
/+(C) = /+(£> = 0,C) = e^f^^ (C) +ib+{()on in terms of the parameters ^T; (; = 1 , . . . , 4), 
/; =/(£, = 0,C = Ki) and the angular velocities J7i = = i?'^), ^2 = l^^^' = J?*'*': 



where 



and 



F-=(-l)"l ,2 



2ini"){(-K„) 



1 -ib+{0 
ib+iO /+(C)/+(C) 



/« -1 

fn ~fn 



(21) 



(22) 



(23) 
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Obviously, the sum of the off-diagonal elements of JV has to vanish, .yVyi + = 0, whence 



1 

1 



4 

n 

«=i 



1+ 



2iJ7(«)(C-^:„) 



(24) 



Since this equation holds identically in C,, one obtains four constraints among f?i , \ ■ 
K4, K2-K3,Ki-K2;fi,... U. Particularly, the 1 /C-term yields 



f[-f\ 

/a ~ It, 



(25) 



Hence the three remaining constraints enable us to express the three similarity variables 
ni{K2-K4), ni{Ki-K2), Qx{K2-Kj,), or, alternatively, ^?2(^3-^4), ^^2(^1 - ^2), 
^2(^2 — Kj,), in terms of the four values of the Ernst potential /i , . . . ,/4. It turns out that 
the constancy of the surface gravity, or, alternatively, the condition A: = on ■e/*, s^^ , gives 
rise to one further constraint. Since the axis values f^{C) determine the solution of the Ernst 
equation uniquely, one may expect a three (real) parameter solution of the two-horizon prob- 
lem (see remark on Eq. dSOll). 



2.4 The double-Kerr-NUT solution 

According to | |2U and i22i . the axis potential /^(C) is a quotient of two polynomials in 
To determine the degree of the polynomials we compare the polynomial structure of the 
matrix elements ( 122b . 

g2f/+(C) = (C-gl)(C-g2)(C-j^3)(C-j^4) ^26) 
^4(0 

i^+(C)=i^, (27) 
/+(C)7^=^, (28) 

/'4(C) 

where p^, 114 are real normalized polynomials of the fourth degree (the fourth power co- 
efficient is equal to one) and the real polynomial P2(C) is of second degree due to ( 125b . 
Replacing /+ and /+ in the third equation by the combination /+ = e^^ -|- ib^ and its 
complex conjugate from the first and the second equation, we get the condition 

[(C-^:i)(C-^:2)(C-^3)(C-^4) + iP2]x 

x[{C-Ki){C-K2){(;-K3){C-K4)- ip2] = -^AiOPA (C) . (29) 

Identifying the zeros of both sides, we see that each bracket of the left hand side has to have 
two zeros of -k^ as well as of p4 (note that the brackets must be complex conjugate to each 
other). Therefore, f^{C) has to be the quotient of two polynomials of second degree, 

(30) 



where the numerator polynomial and the denominator polynomial have the structure 

«2(C)=C' + K + «, 

d2{0=C^ + e( + d, 



(31) 
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and a, b, d, e are complex constants. 

Preparing the continuation of /^(C) off the axis of symmetry into the g-C^ plane (see 
Fig. [T]l we replace these constants by the appropriate parameters d2{Ki) / d2{Ki) and 
h2{Ki)/n2{K-:} , 



Kf + eKi + d 



Pi 



Kf + bKi + a 



' Kf + eKi + d' Kf + bKi+a' 

Equation i |26l l implies c^^^C^i) = and therefore 



a/a; = 1, 



/3/A-=l. 



whence 



f+{Ki) = -r{Ki), i=i,. 



Pi 



In a next step we solve the linear algebraic system 

eK, +d — eaiKi — dcti = Kf[ai — 1), 



/=!,.. .,4, 



(32) 



(33) 



(34) 



(35) 



to obtain e, d and finally rf2(C) in terms of a/, Ki (i = I,. . . ,4). Because of I l34t n2(C) can 
simply be read off from d2{() by replacing a,- by —a; (/ = 1, . . . ,4). Thus we arrive at the 
determinant representation 



/+(C) = 



1 


Kf 






K^ 


Kl 


1 


aiKiiC- 




a2K2{C - K2) 


a^K^iC-K^) 


a4K4{(-K4) 





Ki 




K2 


K3 


K4 





ai(C- 


Ki) 






MC-K4) 





1 




1 


1 


1 


1 








Ki 


Ki 


-1 






OL2K2(C,-K2) 




a4K4{(-K4) 





Ki 




K2 


K3 


K4 





ai(C- 


Ki) 


«2(C-^2) 


a^iC-Ki) 


MC-K4} 





1 




1 


1 


1 



(36) 



for the Ernst potential f{() on the axis ^2/+. 

We will now construct f{g, (). It can be shown that the axis values on .e/+ determine the 
Ernst potential f{g, () everywhere in the g-(^ plane. Hence, if we find a continuation f{g, (^) 
of /^(C) for all £( > and can prove that it satisfies the Ernst equation il2i we have achieved 
our goal. Introducing the "distances" r,- from the points g = 0, = Ki by 



{C-Ki)^ + g^>0, 



,4, 



with the property 



(37) 



(38) 
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, 4) in i36i by r, we arrive at 



Kl 



and replacing the expressions C — Kj{i=\,. 

1 

1 oi\K\r\ aiKiri a},Kj,r}, oi^Kurn 
^1 K2 K4 
a\ri Q-iri OL^r^ ciA^'A 
1111 





K'2 


Kl 


Kl 


a\K\r\ 


a2K2r2 






Kl 


K2 


K3 


K4 


am 








1 


1 


1 


1 



(39) 



A straightforward calculation shows that /, as defined in ( |39b , is indeed a solution of the 
Ernst equation Gad. As we have already mentioned, the remaining gravitational potentials 
k, a (e^^ = can be calculated from / via line integrals. This solution of the vacuum 
Einstein equations represented by the Ernst potential f{g,C) with the axis values OOl l. ( 136b 
is known as the double-Kerr-NUT solution. It depends on seven real parameters: four real 
arguments of a,-, a,ti/ = 1 (i = 1, ... ,4) plus three differences Ki — K2, — K4 ("length" 
of horizons), K2 — ("distance" between the horizons). (Note that the configuration as 
sketched in Fig.[T]can be translated along the (^-axis.) Hence, the solution of the two-horizon 
problem is a (particular) double-Kerr-NUT solution. 



2.5 The equilibrium conditions 

The double-Kerr-NUT solution in the form i |39l > was presented and discussed in [15 | as a 
particular (A' = 2) case of the A'-soliton solution |21 ,220 of the Ernst equation generated 
by the application of N Backlund transformations to an arbitrary seed solution. Applying 
the boundary conditions ^ to the representation ( |39b , Tomimatsu and Kihara I26II de- 
rived a complete set of algebraic equilibrium conditions on the axis of symmetry between 
the parameters q,, Ki (i = 1, . . . ,4). Particular solutions of the algebraic system involving 
numerical results were discussed by Hoenselaers [131, who came to conjecture that the 
double-Kerr-NUT solution cannot describe equilibrium between two aligned rotating black 
holes with positive Komar masses. Hoenselaers and Dietz |8, 12] and Krenzer 1 17] were able 
to prove this conjecture for symmetric configurations (Ki — K2 = K^— K4, f2i = Q2)- 

The final explicit solution of the Tomimatsu-Kihara equilibrium conditions was found 
by Manko, Ruiz and Sanabria-Gomez 1181 . Following their idea, we start with the condition 
/t = Oon i/", 

^±,^2/": k = Q (40) 
and apply it to k calculated from /, see, e.g. |16| . The only condition is 

Qia2 -|- 030^4 = 0. (41) 

' The procedure may seem rather tricky. In fact it reflects steps of the inverse scattering method whose 
explanation is outside the scope of this paper. 

^ A misprint in Ref. 1211 was corrected at the end of Ref. t22l . 
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Combining this result with the two conditions derived from a = on iz/*, iz/" (a again 
calculated from /, see, e.g. |16|) one obtains 



(1+02)^ , (1+ai)^ , ^^23^^24 

-7 ~ 



(42) 



a2 Qi ' Kio,Ki4' 

where 

Kij:=Ki-Kj, f,7=l,...,4. (43) 
Introducing the relative horizon "lengths" 

h = §^, /2 = fi (44) 
A23 A23 

we may express 7, 7' by the scaled quantities Zi , I2 alone, 

_ (l+/2)(l+/l+/2) I+/2 



(45) 



I+/1 ' (l+/i)(l+/i+/2)' 

Setting 

0:304 = —Q1C12 = {aa = I) (46) 

to satisfy ( 141b one obtains the q, (;' = 1 , . . . , 4) from i42[ in terms of the three real parameters 
7, 7' (or, alternatively, /i, h) and arga = </)(« = e"^) 1 18| 



w'a^ +iea a^ + iw'ea 



ai = — : , a2 = - — 

w — lea 1 — iw'ea 

wa^ — a 

Q3 = , 0:4 ■■ 



o? — wa 



1 — wa ' 
where 



(47) 



7'|e(0,l], w:=|v^|e[l,-), £ = ±1. (48) 

Here /j, I2 are arbitrary positive constants and a is a periodic function of a = e"^. 

With the aid of the relations ( 147 l i Manko and Ruiz 1 19| were able to calculate the Komar 
masses Mi, M2 belonging to the horizons J^, J^, respectively, and show that positive 
Komar masses are incompatible with the equilibrium conditions. 

A concise reformulation of the double-Kerr-NUT solution ( 1391 was derived by Yamazaki 



Rn-lRi4-l 

/?23- 1 ^34-1 



air; — air, 
^'7 



7?i2 + l ^14+1 
i?23 + l ^34+1 

whereby the a,-, (/ = 1, ... ,4) have to be taken from ( 147 1 |. Obviously, one can introduce 
dimensionless coordinates q, ( via 

c = ^ (50) 
K23 K23 

and see directly that the Ernst potential, as a function of g and depends only on the three 
parameters li, I2, (f>. We will make use of the formulation ( |49b in the subsequent sections. 
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3 Thermodynamics of the two-horizon solution 

3.1 Thermodynamic quantities 

The best way to get a systematic survey of the relevant physical parameters (state variables) 
of a two-black-hole system and relations among them is to resort to the framework of black 
hole thermodynamics. This theory tells us that the total mass M of the system is a thermo- 
dynamic potential expressed in terms of the independent extensive quantities: horizon areas 
Ai, A2 and angular momenta Ji, J2 of the two black holes. As a consequence of the Gibbs 
formula (see Eq. | |60| | below), the intensive state variables angular velocities fii, and 
surface gravities k\, k2 are functions of the independent quantities too. Furthermore, the 
individual Komar masses M\ , A/2 could play a role. It turns out that all quantities can be 
calculated from the Ernst potential and its derivatives in the points of intersection of horizon 
and symmetry axis (fi = 0, C = AT,-, i = 1, . . . ,4). 

By integrating parts of the Einstein equations over the two horizons Mi , Mi we obtain 
the following relations, 

kxAx=1'k(Kx-^i). k2A2 = 2-k{K^-K^), (51) 

^?lMl = ^(/l-/2), «2M2 = ^(/3-/4), (52) 

niJi = ^--{Ki-K2), f22J2 = ^--iK3-K4), (53) 

where / = f{g = 0, C = K,). Starting with the properties of the Killing vector L = ^ + fi-q 
on the horizon M (see Sec. 12. U one can show that 

Ki-FiJ7i = K2+i^22 = ^/,clc=/f3' (54) 

where /+ and /" are the axis potentials on s^^ and jz/'', respectively. A direct consequence 
of | |5U and l |53l > are the Smarr formulae 1251 

Mi = IQiJi + ^Ai , /• = 1 , 2. (55) 

In order to calculate the ADM mass M, one can make use of the asymptotic behaviour 

/ = 1 - — for r ^ oc, (56) 
r 

where = + (^. Evaluation on leads to 

M=i^lim(l-/+)C. (57) 

Interestingly, the explicit calculation shows that 

M = Mi+M2 (58) 

holds for the 3-parameter solution, i.e. possibly present space-time singularities (see Sec. |5.2| | 
do not contribute to the ADM mass M. As a consequence, we obtain the Smarr formula 



M=Y^^[2niJi + ^Ai) (59) 



for the total mass M. 
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3.2 Gibbs formula 

A regular axisymmetric and stationary vacuum spacetime with black holes obeys the Gibbs 
formula ("first law of black hole thermodynamics") 

<5M = ^(aW, + £-5A,), (60) 

see (U. Hence, the total ADM mass M of the spacetime, as a function of the extensive 
quantities 7, and A, , 

M = M{JiAi), (61) 

is a thermodynamic potential and infinitesimal mass changes 5M between neighbouring 
solutions are given by i l60b . 

There could be spacetime singularities outside the event horizon (see Sec. |5.2| l. Hence 
it is not clear a priori whether the Gibbs formula | |60I > also holds for the special double- 
Kerr-NUT solution ( |49l >, i l47b . Therefore one has to test the validity of ( 160b from the outset. 
For that purpose, we use the formulae from the previous subsection and the Ernst potential 
( |49l >, ( 147 l l to obtain expressions for M, J\,J2, i^i, f^i, ki, k2,A\ and A2. Obviously, all these 
quantities can be written in terms of the four parameters 

{Px,...,P^) = {K2i = K2-K^,w,w\a). (62) 

As an example, the total mass has the explicit form 

M=-^fl + ^). ^ (63) 

2 V w / M + esm^cosi;/) 

with 

M:=\ + -\w' + — \sm4) )cos<^. (64) 

2\ w' J 2 \ w J 

Equation ( I6OI I is equivalent to the four equations 

_ = + + + ,= i,...,4. (65) 



A straightforward calculation shows that ( I65l > is indeed satisfied. Therefore, we may con- 
clude that possible singularities do not contribute to &M and the first law of thermodynamics 
(Ho} holds. 



4 The sub-extremality of black holes 

Following Booth and Fairhurst 1 6 1, we will assume that a physically reasonable non-degener- 
at^^ black hole should be sub-extremal, i.e. characterized through the existence of trapped 
surfaces (surfaces with a negative expansion of outgoing null geodesies) in every sufficiently 
small interior neighbourhood of the event horizon. It can be shown |9 1 that any such axisym- 
metric and stationary sub-extremal black hole satisfies the inequalit)0 

87r|7|<A, (66) 

' The degenerate (extremal) case requires special attention. 

^ Note that the inequality j66t can be generalized to the Einstein-Maxwell case, i.e. to electrically charged 
black holes, see 1101 . 
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i.e. for given event horizon area A, there exists an upper bound for the absolute value of the 

angular momentum |/|. 

In order to test explicitly whether the two gravitational sources with the horizons Jfj , 
(tentative black holes) in the double-Kerr-NUT solution l l49t . ( 147 1 satisfy this inequality, 

we calculate the quantities 

P,:=^, (67) 

"■I 

We obtain the remarkably simple expressions 

l+<?>w' w(w — 

Pi=e-— - — rr, Pi = e— — , (68) 

w'(^' + w') 1— 

where 

cos (?!) + £ sin e [-^2,^2]. (69) 
Hence, the inequality 1 I66I 1 for each of the two black holes is equivalent to 



w'2 + 2$w' + 1 



and 



respectively. Taking into account the allowed parameter ranges w 6 [1,°°), w' 6 (0, 1], these 
inequalities can only hold if 

w'^ + 2#w'+l<0 and w^-2#w+l<0. (72) 

However, this implies $w' < and <Pw > in contradiction to w' > and w > 0. 

Thus we have proved the non-existence of a stationary and axisymmetric two-black- 
hole configuration with separate horizons (see Fig. [T]l, i.e. the spin-spin repulsion of two 
aligned black holes cannot compensate for their gravitational attraction. The non-existence 
theorem is essentially based on the inequality ( I66I 1 that is as shown in ^9] a consequence 
of a defining geometrical black hole property (6]| (the case of extremal black holes requires 
special attention). 



5 Further properties of the double-Kerr-NUT solution 

As we have seen in the previous section, the equilibrium of two aligned black holes is im- 
possible. The only candidate for a solution of the balance problem — the double-Kerr-NUT 
solution — has to be dismissed as a physically irrelevant solution as discussed above. Never- 
theless, it is interesting to study further properties of the solution l |49l >, i |47| |. In the following 
two subsections we comment shortly on the interior of the two gravitational sources and 
give numerical evidence for the existence of singularities outside the horizons and J^. 
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5.1 The interior of black holes 

It was shown in 1 2 1 that every axisymmetric and stationary black hole, which is regular in an 
exterior neighbourhood of the event horizon, also possesses a regular interior region inside 
the event horizon. In particular, there always exists a regular inner Cauchy horizon and the 
inner solution does not develop singularities before this horizon is reached. Moreover, the 
spacetime is even regular at the Cauchy horizon, provided that the angular momentum / of 
the black hole does not vanish. Remarkably, the areas A and A'^^ of event and inner Cauchy 
horizon satisfy the equatior0 

{8TTjf = A™A. (73) 

It is interesting to test this relation explicitly for both of the two gravitational sources in 
the double-Kerr-NUT solution i4% . i47i . For that purpose, we calculate the areas M , J^i 
of the horizons J^i, J^2 using iSli and i |54| |. As shown in (2], the areas of Cauchy horizons 
can be calculated from analytical continuations of and into regions with C, < Ki and 
< K-}, respectively. In this way we obtain 

Using these formulae, the explicit calculation shows that the equations 

{S,TvJif =Af^Ai, (=1,2, (75) 
are indeed satisfied, i.e. ( 173b holds for both gravitational sources. 



5.2 Singularities outside the black holes 

As we have proved in Sec. |4l at least one of the two "black holes" in the double-Kerr- 
NUT solution i |49b , ( 147b is not sub-extremal for which reason the solution is not physically 
reasonable. It may then well be that singularities outside the horizons appear 

To tackle this problem we ask whether the determinant in the denominator of the rep- 
resentation ( |49b of the Ernst potential / has zeros off the axis. A numerical study for a 
large number of parameter values shows that this is indeed the case, see Figs.|2]and[3] As a 
consequence, the Ernst potential becomes singular at these zeros (it can be shown that the 
numerator does not vanish at the same coordinate positions), i.e. there exist singular rings 
outside the horizons. Our numerical investigations seem to indicate that every double-Kerr- 
NUT solution with the special parameter relations ( 147 b suffers from the presence of singular 
rings. 



6 Summary 

The stationary equilibrium of two aligned rotating black holes can be described by a bound- 
ary value problem for two separate (Killing-) horizons (see Fig. [T]l. Applying the inverse 
(scattering) method, one can show that the solution of the problem is a (particular case 
of the) double-Kerr-NUT solution (a solution originally generated by a two-fold Backlund 

^ Note that these statements can be generalized to black hole spacetimes with electromagnetic fields, see 

i3Tm . 
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Fig. 2 Singularities of the Ernst potential: The plots show, for two different configurations, curves along 
which the real part (solid curve) and imaginary part (dashed curve) of the determinant in the denominator of 
the Ernst potential / vanish. At the intersection points Si, S2, (53 ), the Ernst potential diverges. The horizons 
Jfl, ,S^2 are marked as black lines on the f-axis. 

Parameters: = ^tt, w = 10/3, w' = 0.3 (left panel) and </> = —0.1, w = 1.3, w' = 0.5 (light panel). 




Q 

Fig. 3 Singular Ernst-Potential: Real and imaginary part of the Ernst potential / for the example configuration 
in Fig.|2] right panel. 



transformation of Minkowski space). The regularity conditions to be satisfied by the metric 
on the axis of symmetry outside the two horizons restrict the number of free parameters 
entering the solution. The resulting 3-parameter solution (written in dimensionless coordi- 
nates) does not satisfy the characteristic condition 87r|7| < A for each of the two black holes 
(/: angular momentum, A: area of the horizon). Since this inequality is a consequence of 
the geometry of trapped surfaces in the interior vicinity of the event horizon of every sub- 
extremal black hole, there exists no stationary equilibrium configuration for two aligned 
sub-extremal black holes. 
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